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Abstract. Let if be a totally real Galois number field. C. J. Hillar proved 
that if / G Q[xi, • • ■ , x n ] is a sum of m squares in K[xi, ■ ■ ■ , x n ], then / is a 
sum of 2V(m) squares inQ[a;i, ...,x n ], where N(m) < 2[ K:< M +1 • ( [X: ® 1+1 ) -4m, 
the proof being constructive. 

We show in fact that N(m) < (4[K : Q] — 3)-m, the proof being constructive 
as well. 

1. Introduction 

In the theory of semidcfinite linear programming, there is a question by Sturmfels 

Question 1.1 (Sturmfels). If / S Q[xi, . . . , x n ] is a sum of squares in R[xi, . . . , x n ], 
then is / also a sum of squares in Q[xi, . . . , x n ] ? 

Hillar ([3]) answers the question in the case where the sum of squares has coef- 
ficients in a totally real Galois number field : 

Theorem 1.2 (Hillar). Let f £ Q[x±, . . . , x n ] be a sum ofm squares in K[x\, . . . , x n ] 
where K is a totally real Galois extension ofQ. Then, f is a sum of 

IK . oi + i f[K:Q] + V 
2 



4m 



squares m HJ>[£i, . . . , x n \. 

The aim of this note is to show, modifying a little bit Hillar's proof, that only 
(4[K : Q] — 3) • m squares are needed (that is Theorem 13. L[) . Moreover, as in [3], 
the argument is constructive. 

2. Hillar's method 

Having in mind the Lagrange's four squares Theorem, we focus ourselves on 
rational sum of squares i.e. linear combination of squares with positive rationnal 
coefficients. 

Let if be a totally real Galois extension of Q which we write K = Q(0) with 9 
a real algebraic number, all of whose conjugates are also real. We set r = [K ; Q] 
and G = Gal(if/Q). 

Let / € Q[xi, . . . , x n ] be a sum of m squares in K[x\, . . . , x n ], namely / = 
SfcLi /fei with fk € K[x\, . . . ,x n ]. Summing over all actions of G (i.e. "averag- 
ing"), we get 
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1 m 

Next, we write each fk in the form 



G 

1 k=l <J£G 



r-1 

fk = J2 

i=0 

where qi € Q\xx, . . . , x n ]. Then, 

r /r-1 \ 2 

(2) EH) 2 = E E*M)' 

(TEG 3=1 Vi=0 / 

We may write this sum of squares as the following product of matrices 

' ' I <j x ... {axOy- 1 \ ( q 
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\ 1 a r d . 


.. (^ey- 1 




1 a r 9 ... {a r 9y- L / \ g r _i 



We obtain what is called a Gram matrix (cf pQ) associated to the sum of squares 
in ©. Let 



1 <7i0 ... {<J X ey- 1 \ T ( 1 <7i0 ... (o-i0) r 



G = 



i a r 9 ... (drey- 1 J \ i o>0 ... {a r ey- x 

Note that the entries of G are in Q since they are invariant under the <Tj's. 

Now, we come to the slight modification of the proof of Hillar that will improve 
the bound. 

3. LU-DECOMPOSITION OF THE GRAM MATRIX 

If u{x) denotes the minimal polynomial of the Galois extension K over (Q, then 
the (i, j)-th entry of the matrix G is the i + j — 2-th Newton sum of (oi, . . . , o>) 
the roots of u(x). It is well known that the rank of G is equal to r and its signature 
(the difference between the positive eigenvalues and the negative ones) is equal to 
the number of real roots of u(x) (see for instance [2 Theorem 4.57]). In our case, 
we readily deduce that G is a positive definite matrix since K is totally real. Thus, 
all its principal minors are different from zero (they are strictly positive !) and G 
admits a LU-decomposition which we may put in a symmetric form 

G = U T DU 

where D is diagonal and U is upper triangular with diagonal identity, and U, D 
have rational entries. 

We may view this decomposition as a matricial realization of the Gauss algorithm 
which reduce the quadratic form given by G. 
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Now, if we denote by f±, . . . , f r the polynomials in Q[x%, . . . , x r ] appearing as 

( »\ 

the rows of the matrix U x and by d\ , . . . , d r the rational entries onto 

\ Qr-1 J 

the diagonal of D, then we get from ([2]) the identity : 

(3) X>/*) 2 - ytfl + ■ ■ ■ + ^-/r 

This construction leads to 

Theorem 3.1. Let f e Q[xi, . . . , x n ] be a sum ofm squares in K[x\, . . . , x n ], where 
K is a totally real Galois extension of Q. Then, f is a sum of {A[K : Q] — 3) • m 
squares in Q[xi, . . . ,x n ]. 

Proof. By ^ and (|3|), it suffices to apply Lagrange's four squares Theorem to get 
that / is a sum of i[K : Q] ■ m squares in Q[xi, . . . , x n \. 

But let us note that the first diagonal entry of D is always d\ = r = [K : Q] . 
Then, by the averaging process the first coefficient appearing in the rational sum of 
squares in (J3J) is -^j = 1 : already a square in Q ! Whereas the others coefficients 

t^It in the rational sum of squares could be any positive rational which we rewrite 
as a sum of 4 squares. This concludes the proof. □ 

Remark 3.2. Beware that if we perform the Cholesky algorithm to the matrix G 
instead of the LU-decomposition, it yields a factorisation G = U T U where U is 
lower triangular but with entries in Q[\/di, . . . , for some integers d\, . . . , d r . 
Then, an averaging argument would produce identities over Q but will raise the 
number of squares by an unexpected multiplicative factor 2[ A_:< ®. 

Let us consider as an example, the simple case of quadratic extensions : 



Example 3.3. Let K = Q(vu) where d £ Q is not a square. The extension K is 
always Galois, and it is totally real if d > 0. 

Let / S Q(xi,...,x n ) be such that / = J2k=i( a k + hVd) 2 with ak,b k e 
Q(xi, . . . , x n )- Since / has rational coefficients, by averaging we get 

-. m m 
k=l fe=l 

It remains to write d as a sum of I < 4 squares of rationals, and we get that / is a 
sum of at most (1 + I) ■ m squares in Q(ii, . . . , x n ). 

As another illustration, we apply our method to [H Example 1.7] : 

Example 3.4. Consider the polynomial 

/ = 3 - 12y - 6x 3 + 18y 2 + 3x 6 + I2x 3 y - 6xy 3 + 6x 2 y 4 

which is the following sum of squares 

/ = (x 3 + a 2 y + [3xy 2 - l) 2 + (x 3 + /3 2 y + jxy 2 - l) 2 + (x 3 + j 2 y + axy 2 - l) 2 

in Q(a)[x, y] where a, /?, 7 are the real roots of the polynomial u(x) = x 3 — 3x + 1. 
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We do not need to average and directly compute the matrix G and its symmetric 
LU-decompositon 

3 6 \ / 1 00\/300\/l0 2 \ 
6 -3 = 10 06 01 -| 
6 -3 18 / V 2 -\ 1 / V I / V 1 / 

Because of the relations (3 = 2 — a — a 2 and 7 = a 2 — 2, the vector of polynomials 
<Z = (<Zo, Si) Q2) T is q = {x 3 + 2xy 2 — 1, —xy 2 , y — xy 2 ) T and hence 

/ = 3 (V + 2xy 2 - 1) + 2 (y - xy 2 )) +6 (-xy 2 - -(y - xy 2 )) 

a rationnal sum of 3 squares, to compare with the rational sum of 6 squares 
obtained in [3j. 
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